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ABSTRACT 


The problem of aerodynamic interference between axisymmetric 
bodies in inviscid, incompressible flow has been studied by 
separating the flow into axial flow and lateral flow. The 
problem of axial flow has been solved by representing isolated 
bodies in free stream by continuous source distribution of 
polynomial form along the axes of the bodies and the interfering 
flows by an additional continuous dipole distribution of 
polynomial form after Milne -Thompson's circle image theorem. 

The solution for the lateral flow is obtained by invoking slender 
body theory and using the method of two-dimensional multipole singu 
larities (MPS ) . Flow tangency conditions over the surfaces of 
the body has been used for the axial flow problem, while 
constancy of stream function on the surfaces of the bodies has 
been used for the lateral flow. This method has been applied to 
study the interference between a pair of parabolic bodies 
of various thicknesses and for varying proximities^ • The per- 
turbation velocities, the pressure distribution and the side force 
co-efficients and moments are computed and .presented . This method 
of solving interference problem by using dipole distribution 
in axial flow and two-dimensional MF5 in lateral flow appears to 
be a rapid and fairly accurate method for preliminary calculation 
of aerodynamic loads. This facilitates more efficient use of 
panel methods if required for greater accuracy. 
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CHAPTER I 


Introduction 

i r> 

This study is about aerodynamic b ehaviou r of external 
stores and armaments carried by military aircraft . The 
interference could be between wing and the external stores, 
fuselage and external stores, or between the external 
stores themselves. The problem of interference between 
external stores, which are mostly axially symmetric bodies, 
has been taken up for the present study. The solution of 
this problem has been obtained by s eparating the ' axial flow' 
and the 'lateral flow'. The 'axial flow' solution is based 
on method of images and by a continuous source distribution , 
represented by a polynomial function, displaced from the 
axial line towards the interfering' body. Lateral flow has 
been solved using slender body assumptions, treating the 
flow at various cross-sections along the axis as 2-dimensional. 
Multipole singularities like doublets quadrupoles etc. 
have been used for nullifying the cross-flow. 

1 * 1 Motivation for the study 

Modern military aircraft carry large amount of 
external stores and armaments like bombs, missiles, rockets 
and drop tanks under the wing and under the fuselage. Most 
of these stores are carried in close proximities to each 
other. There is considerable mutual interference of flow 
between these bodies. Proper knowledge of interference 
is necessary to calculate aerodynamic forces and moments 



2 . 


for safe carriage, launching and accurate delivery of these 
stores to their targets. An approximate, yet a rapid method 
of estimation of these effects is helpful in optimising 
and utilizing the available under the wing/fuselage area 
for maximum armament carrying capacity. 

1 . 2 Literature Survey 

The subject of flows past system of axi-symmetric bodies 
has been studied by many authors. Most of them have used 
singularities like sources, sinks, doublets and still 
higher derivatives of these singularities called multipoles 
to represent various shapes and bodies. These singularities 
could be distributed within the- body or on the surface of the 
body. The former is called ’Internal Singularity’ method 
and the latter 'Surface Singularity' method. The ’surface 
singularity' method is also related to the 'Panel method' 
wherein the surface is discretised into large number of 
panels. 'Panel method* is accurate but its accuracy 
depends on number of panels. So for greater accuracy, larger 
number of panels and hence larger computations are required. 
The method of internal singularity is simpler, rapid and 
reasonably accurate for smooth bodies. 

The pioneering work in flows over axi-symmetric bodies 
was done by Rankine [l] , who represented the axisymmetric 
body by source-sink pairs, along the axis of the body. 
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The bodies which may be represented by source-sink 
combinations are called 'Rankine bodies' or 'Rankine ovals' 
Von-Karman [ 2 ] represented a flow past body of revolution 
by continuous distribution of sources/sinks along the axis 
of the body. He used different source strength for 
different segments but the strength being constant over 
a particular element. For lateral flow he used continuous 
distribution of doublets, the doublets aligned opposite to 
the undisturbed lateral flow. James [3j obtained an 
analytic solution of velocity potential 0 for flow over an 
axi-symmetric body by conformal mapping of the meridian 
profile of the body. He transformed meridian geometry of 
body into a circle and then approximated the velocity 

♦ 

potential in terms of elementary functions such as Fourier 
Series, Legendre or Chebyshev polynomials. He obtained 
results which are comparable to those obtained by 'panel 
method . ' 

Zedan and Dalton [4] represented the body by parabolic 
variation of source strengths along axial segments, of the 
body. By using 12 elements of source strengths each of 
2 or 3 degree polynomial, they obtained accurate results 
for representing even thick bodies of fineness ratio upto 
2, by iterative method. Their method is suitable for blunt 
nosed bodies also. Campbell [5] obtained least square 
solution of discrete distribution of sources along the axis 
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His method is unsuitable for blunt bodies and also the number 
of sources required are large. Christopher and Shaw [6] 
used multi-polynomials of degrees 8 to 10 to represent bodies 
They obtained accurate results for thick as well as blunt 
bodies . 

Martin, Smith and Saunders [7] studied interference 
between two identical a xi- symmetric bodies by using discrete 
distribution of sources for representing original bodies 
and using discrete distribution of dipoles along a displaced 
line which may be called imago or delta line- for the 
interfering flew. The drawback with discrete distribution 
of sources is that boundary conditions are satisfied at 
control points but there is a large error at points other 
than at control points, especially near nose or tail of the 
body. Christopher and Shaw [8j represented the bodies by 
multi-polynomials and used 3-dimensional multipole singular! 
ties (iviPS) for the cross-flow due to the interference. 

They used MPS i.e. doublets, quadrupoles and octopoles 
each represented by polynomials of degree 8 to 10. This 
method has two drawbacks; firstly it has large number of 
unknown co-efficients (as large as 120 ) because of multi- 


polynomials and MPS. 
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1.3 Outline of the present work 


In chapter 2, the axial flow problem is formulated. 
Continuous source distribution represented by 8-degree 
polynomial is used for isolated bodies. For interference 
an additional continuous distribution of dipole strength 
function represented by a 8 degree polynomial is used. 

The mathematical Expressions of velocity potential 0 
and perturbation velocity components i.e. v^, u x , are 

given along with governing equation and boundary conditions. 
The axial flow over a pair of axially- symmetric bodies, 
laterally separated, are studied for varying thickness ratios 
and proximities. Results are presented. In chapter 3, 
the lateral flow problem is formulated using slender body 
theory. The concept and application of 2-dimensional MES, 
is explained. The mathematical expressions of stream function 
and perturbation velocity components Vq^ and are 
given. Boundary conditions are also described. The 
two types of lateral flow viz., yawing and incident are 
solved separately and results are presented. The aerodynamic 
load distribution Cp, forces and moments are obtained. 


In Chapter 4, the results are discussed. Limitations 

of 

of discrete distribution/sources and factors responsible 
for err->r are stated. Suggestions for further extension 
and utilization of present work are given. 



CHAPTER 2 


6 . 


Axial Flow 

2 • 1 Potential f low' solution 

This chapter is devoted to analytic work and the 
mathematical formulations necessary to obtain solution of 
velocity potential 0, and to determine polynomial functions 
f (x)if ( x) of the source and- dipole distributions, the per- 
turbation velocity components (u , v , v~ ) , pressure 

X X' V Ci 

coefficient Cp, normal/side force co-efficients ^z, Cy and 

the moment coefficients C . C . 

mz my 

. Subsonic flows, at low mach No. past axi-symmetric 
slender bodies may be treated as incompressible and where 
Reynolds he. is near 2x10° or greater (Turbulent Boundary 
Layer) may bo apprximatec to inviscid flows [ 6 j because the 
boundary layer is thin and not subject to separation. Therefore, 
axi-symmetric bodies at small angles of yaw or incidence 
(limited to 6°) may also be treated by this approach. Solution 
for such flows have been obtained by several authors by 
treating them as inviscid incompressible flows. The 
governing equation for the velocity potential 0 for such 
flows is the weli-kncwp ’Laplace Equation 1 . 

? 2 0 = C 

0 +- 0+0 =0 
^xx r *r *rr 


2.1 
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The s.lutiun to the above Governing equation should £ 
further satisfy the boundary conditions as given here. 

1. 'Flow Tangoncy' Conditions V , n = 0 on the body (2.2) 
The flew on body surface is tangential and the velocity 
normal to surface is zero. Since V = V0 so V0.n=O. 


2. 'Remote Flow' conditions s- Perturbation velocity potential 
at infinity is zero 

0 = 0 at r = co (2.3) 

An axially symmetric body is represented by r = r(x) 
r - r(x) = f(x,r) = 0 

n = Vf = c^Eix}_ i + j 

where i,j are unit vectors in x and r directions. 

So V. n J = vZ.n :< = /U M + u x ^ x - v r = 0 (2.4) 


Where u and v are the perturbation velocity components 
x r . 

in x anc! r directions. 


From 2.4 we can state 


Total radial vel. 
Total axial vel. 


= slope- of body. 


(2.4) (b) 


2 . 2 Isolated body solution 

As indicated earlier in chapter 1, an isolated axi- 
symmetric body in uniform flow is represented by internal 
singularities distributed along the axis of the body. 
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The perturbation velocity potential 0 = f ~ (2.5) 

o fR 

Vvhere f(x) is the source distribution function along the 

axis of body. 

without loss of generality the body length is taken as 1 
R = (x-x ) ^ + r ^ , the distance P D Ref. fig.(l) 

P = Field point P (x,r»0) 

D = the position (x,oj of linear element of singularity 
dis tribution 

x = dummy variable, which represents location of 
singularity distribution. 

The axial and radial perturbation velocity components at 


a field point P are obtained by differentiating 0 along 
respective directions. 


v r 

II 

r . 

= 

■I 1 

0 

f(x) dx 



(2.6) 

u 

X 

_ 

dx 

II 

4^1 1 

cnjM 

. f 1 

%/ 

0 

f ( x ) (x - 
R 3 ^ 

X ) dx 

>9 

(2.7) 



= X- 

I 1 

0 

f(x) dx 
r 3/2 

- / 1 ' 
0 

f ( x) x dx 
-^372” 

(2.8) 

2 • 3 

Solution 

of polynomial source 

distribution f(x) 



In order to represent the body accurately, a review 
of work of previous authors and their choice of f(x) is 
essential. Accuracy and sensitivity of the method is dependent 



9 . 


upon proper choice of f(x). Von Karrnan [ 2 ] used f(x) constant 
over various segments along axis. Zedap [3] used linear 
variation of f(x) and in a later work used parabolic variation 
of f (x) over different segments. Christopher and Shaw [6] 
chose multiple polynomial, (3-polynomials of degree 10 each). 
Based on their experience, we have tried and chosen a single 
polynomial. A single polynomial of degree 8 is found to 
represent an axi-symmetric body fairly accurately 

So f(x) = C q + ^ x + C 2 x 2 - - - ~ u n x 

n 


= Z 

o 


C . x J 


(2.9) 


Solution for f(xj needs evaluation of the integral of form 

.1 


T r x x^dx 

1 3 = { -7372 


( 2 . 10 ) 


R = R (x, x- , r»$) 

It is interesting to state hc-re that Ij is easily deter- 
mined from the table of integrals and it is not calculated by 
quadrature numerical integration. The expressions for the 
indefinite integrals of the form of I . are given in recursive 
form in the book 'Table of integrals, Series and Products' by 
Sradshtcyn and Ryzhik [9]. These are reproduced in App. A. 

Cj is evaluated by imposing the boundary conditions .at n+1 
control stations P on the body. 


Substituting v 


r 5 u x 
= JU + u v 

X 


in eqn. (2.4) viz. 
dr (x ) 

x dx 
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whe re 


dr(x) 
dx 


= slope of the body at field point P 


Then eqn. (2.4) in vector form is 


n 

j=o 


n n dr(x) 

c : x b = u ~ + x i - =0 c j xI j T + c j x h+i x dT 


Taking = 1, (2.1l) gives (n+1) eqn. by moving field point 
P^ at n+1 control point of the form 
AC = B 
C = A“ 1 B 


where C = unknown co-ef f icients C - 


A 


= Sq. matrix of order (n+l) = A(i,j) = r^xl . 


dr. 

■ (x i xI j - I j+i )x 


dr. 

8(x) = 


By solution of eqn. (2.11) > we obtain (n+l) values of 
unknown variable Cj knowing C j , the source distribution 
function f( x) , and the flow variables and the perturbation 
velocity components v (i), u (i) are readily computed. 

J- }\ 

2 . 4 Flow with interference 

Consider two axi-symmotric bodies b^ and b^ laterally 
separated by a distance ' a' and axially placed in a unifom 
flow U w . (fig. 2) . . 


Due to the presence of body b^? the body b^ is no longer 
placed in a uniform flow hence it may not be adequately 
represented by the source distribution along the axis. Similarly 
body b 2 cannot bo adequately represented by a source distri- 
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bution along the axis. According to circle image theorem of 
kilne-Thompson [10]^ £he source distribution have to be 
displaced off the axes towards each other as shown in fig. 2(a) 
and 2(b) . The lines on which sources are distributed may be 
called as image or delta lines. Iviartin [7] has given an 
expression for calculation of displacement distance d , for 
a 2-dirnensional case, for two identical bodies by making 
repeated applications of circle image theorem. 

6 = f- x -ifl - (^L.) 2 ? (2.12) 

whore 6, r arc the displacement or delta distance, radius 
of circles and C 2 at any cross-section cc. From this 
equation it is clear that 6 depends on geometry and not on 
strength of sources i.e. ^ = 6 (r , a). 

martin [7j has solved the interference problem by 
distributing discrete distribution of dipolos i.e. sources of 
the dipole along delta line and the sinks along the axis of the 
body ref. fig. (2b). The discrete distribution of singularities 
gives pressure distribution which is in good agreement at 
control points but not at test points. Therefore, a continuous 
distribution of dipole f^( X) is tried and utilized to nullify 
the flow induced by the presence of the other body so that 
the resultant flow is tangential to body surface (the source 
distribution of dipole along the delta line and the sink 
distribution along the axis) . The resulting flow is obtainable 
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by superposition of the flow due to source strength functions 
f-^(x), f 2 ( x; representing isolated bodies b^ and b^ uniform 
free stream and that due to dipole functions f^d(x), f^dCx). 
The tangoncy boundary condition given by eqn . 2.4 is applicable 
fur the total flow. i.o. 

dr (x) 

(2.13) 


rt 


Where 


= (U + u . ) 

03 ’ Xu 


X 


GX 


are the total or resulting radial and 


rt' ’ “xt 

axial perturbation velocity components due t..< all singularities 
present in the flow. 


(2.14) 


V rt 

II 

< 

CO 

+ v s2 

r — 1 

> 

V d2 

u xt 

= U , 

si 

+ U s2 

* u ai 

+ U d2 

u hero 

v and u 

represent 

radial 

and axial 


of velocities subscript s,d, represent these components due to 
sources and dipoles. Substituting values of u^ and v 
from eqn. (2.14) in eqn. (2.13) we obtain an )thor relation 
given as below 


A V, 


dr (x) 


Au, 


where *v 


Au 


x 


(2.15) 

s are increments in radial and axial 


induced velocities over and above those '.-f single body in 


uniform flow. 
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It follows fr^m oqn. (2.15) that the induced velocities due tc 
additional singularities i.e. (x)> f-^x), f 2 c’(x) also 
satisfy tangency condition on field point P. , on the body b-^. 


4 v r ‘ 

CM 

CO 

> 

11 

• + v di 

V cl2 

A u x 

CM 

£0 

3 

II 

+ u di 

+ u d2 

Expros 

sions 

for v di’ 

V d 2 ? u dl and u d2 are given in 

. (2.17) 

and 

(2.18) . 



2.5 Solution for the dipole strength function f d (x) 


The velocity potential 0 i duo to a dipole distribution 

U. 

at a field point ,p(x ,r,©)is given as 


^d = + $ 


cip 


dn 


-1 r 1 f d^ x ^ f r 1 f d^ x l ^ 

4 rt r j " wr 


4it 


(2.16) 


where 0 dp , 0 dn are velocity potential due to source and sink 
respectively. 


Her = ( x-x ) 2 + ro- 2 Ref. fig. (3) 


o 


=. 4- 

Let v.,^) = radial velocity due to dipole fn.f,,(x) at P (x,r-,©) 

Ci 


0. 


r 




rc 


4 ?r 




f, ( x) dx 


' Q 

R 


~3j2 


f 


1 f,(x) dx 


4% 


R 


w 


(2.17) 


u d = axial perturbation velocity due to dipole function 


f d (x) at P (x. , r ,©) 




“3^ 


i _ _ i„ r 1 f d ( ^ } (x ~ x - } dx 4. a, r 1 f d 

c 4 It 0 R^ 3 / 2 47t 0 


1 f r1 (x)(x-x ) dx 


■^72 


(2.18) 
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Substituting the value of in equation (2.15) the 

solution for the (n+l) co-efficients of polynomial function 
f ,(x) is obtained by moving fiolcl point P(x ,r,©)at (n+l) 
control points. The exact form of equation (2.15) in vector 
form are given in App. ’ B * alongwith solution. 

2 . 6 Determination of aerodynamic coefficients 


Pressure distribution must b 


finding other co-efficients like 

dc. 


V 


Cp = 1 - 


loading -y— 1 

o. A 

n 


momenx 


U 


■ 2 


e established first before 

side force co-efficients 

coefficient C, Cnd^etc. 
myj 


v = 
n 



u 

v. 


rt 

rb' 



rb 


v. 


rb 


Resultant velocity 

(U + u 2 ’ + v , 2 + V Q 2 ) 

w x rb Q a J 

= sum of perturbation velocity components in x direofcon 
= sum of perturbation velocity components inradial direc 
tion. 

= perturbation velocity components tangential at Field 

Point P obtained by resolving v r due to f^(x), 

fj_d(x), f^dOx). 

d (x) 

= (U^+u^) x according to slender body theory 

and B.C. 

dr ( x) 

= Approx. — x ^oo ^- n case of u . being small, 

which is the case when single slender body is placed in 
uniform stream. But u , cannot bo ignored in case of 

X 

interfering flows in close proximity. 

Knowing Cp other co-efficients are determined as discussed 



below. 


15 . 



dx 



f 2n c P (6 > x) 


cos © d© 


Vvhere S = ref. area = 2 x r 


( 2 . 20 ) 
( 2.21 ) 


clc 


dx 


r 

T 


27V 


— J Cp(©,x) sin © d© 
x o 

Then side force co-efficients 


C y = / 


dc 


X. 


dx 


The normal force coefficient 
dc ' . 

C = / — a:. 1 " dx 

o 


dx 


ivioment co-efficients 
1 dc 




■nd (+x) dx 


C = J 

0 


1 dc 

~sr ( -x) dx 


( 2 . 22 ) 


(2 .23- 


(2.24) 


(2.25) 


All these co-efficients wore determined by numerical integra- 
tion by Simpsons rule dividing periphery in 24 intervals 
and the axial body in 10 intervals. 


2 . 7 Presentation of results 

The proposed method of representing the axially symmetric 
bodies by polynomial distribution of sources and dipoles, is 
applied to the study of flow ever a family of sharp-endec! 
bodies generated by arcs of parabola. This family of bodies is 
chosen in order t_ compare the results of the proposed method 
with the results of other authors who presented their results 


for similar bodies. 



The variation of b as function of x and PR are plotted 
in fig. 9. The source, and dipole strength functions f(x) and 
f ,(x) are plotted in fig. 10. The results of variation in error 
as function of x,6, 17 and proximity (PR) are plotted in 
fig. 11, 12 and 13. The perturbation velocity components u x , 

V~ are plotted in fig. 14,15 and 16. Pressure distribution 

C7 Q, 

(Cp) results are plotted in fig. 17. The side load distributions 
dc 

” — ^ and side force co—ef fieitts C v are plotted in fig. 18,19 and 
dx y 1 * 

20 . Results of C , and C m „ are not plotted as these values are 

mz my 

just half of C z and C y for the family of bodies chosen i.e. 
symmetrical bodies about mid point. 
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CHAPTER 3 
Lateral Flow 


This chapter deals with the application of multipole 
singularities (MPS) for solving the lateral flow. When the 
onset flow is not axial, there will be a lateral component. 

This lateral flow may further be resolved into the yawing and 
incident flows, the yawing flow being the component along the 
line joining the centres of the bodies and the incident flow 
being normal to it. The system of two bodies may then be 
considered to be at an angle of yaw or at an angle of incidence 
or both. (Ref. fig. 5(b) and 6(b)). 

Vve propose to determine the general flow over the axially 
symmetric body by application of slender body theory, treating 
each cross-section as 2-dimensional. This is valid for flows 
at an angle of incidence or yaw not exceeding 6° over slender 
bodies. 

3.1 Interfering flows 

Consider the flow as shown in fig. 5(a) and 5(b). If a single 
body C^ was present in uniform stream V^, then a doublet of 
appropriate strength aligned in y-direction, the direction of 
onset flchw, could have satisfied the boundary conditions. Then 
the solution was v ei = " 2 v - sin w, V r = 0 on body surface. 
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But the- presence of body C r . makes the onset flow V non-uniform. 

2 “flow 

It is proposed to represent the circular body in non-uni form ^ey 

a series of MPS placed at the centre of the circle. The stream 

function of a MPS of order N is obtained by successive 

differentiation N times of stream function of a source along 

a direction. The direction of differentiation may be called the 

direction of alignment of the MPS. The mathematical expressions 

for , and V r induced by a M_PS are given in the next section. 

3.2 . Stream function of multipole singularity 


Consider a uniform onset flow (Ref. fig. 7) 

If a source of strength ivi is located at A (0,0), then 
stream function^ at a field point 
P (r,©) or P (y> z )auo to the source at 



in polar co-ordinates 

M s 

O = 7T” ©. 

Ts 2x 


3.1 


■where ©=tan 



d©_ -z -sin & 

dy“ r 2 r 


r = Y (y 2 +z 2 ) 





denoting y; Q as ''f' due to doublet, quadrupole and octopoles 

aligned in y direction by successive differentiation with 
reference to y direction we get. 


V'D 2ftr 


sin © 


sin 2 © 


= - 2 / 


sin 3 © 


So^n = Stream function due to nth order differentiation. 


we have 


Li/ = (-1) (n--l ) | M ; n sin n © 

T n _ n 

2'C r 


(3.2) 


If a series of MPS are placed at A(o,o) then stream function 
V-'mp due to all those MPS uptoordor N at a field point P (r,©) 
is given by the following relation 


^mp 


' D + fei + fb + 


= E An Jvi n sin n© 

n=l n 

r 


(3.3) 


where An 


2rc ' 


Co-efficient A n can be merged with strength M n . 


2 A 


m(n) 


sin n© 


where A m (n)= A n x fv -n * 


It may be observed from thoabove eqn. (3.3) thatvp^ i.e. 
due to series of nth order MFC, aligned in y direction is sum 
of all sine terms of a Eburier scries in 9. This implies that 
H' mp ^"ti-symmetric about y axis. 
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3 . 3 Yawing flow 

The solution to the problem of lateral flow (yawing or 
incident) over the circles and is attempted by represen- 
ting each circle and C 0 by a series of multipoles of different 
order say, up to N, aligned along the direction of flow (yaw or 
incident) and placed at the centres of circles. The strength 
of 2N MPS are determined by imposing the boundary conditions 
^'=const. at N control points on each circle. An alternate 
way of imposing boundary condition is V r =0 at N control points 
on each circles. 


From the boundary c onditions the total stream function, due 
to onset flow and both the MPS at the centres of circles C 1 and 
Cp, must bo zero. Then resulting equation j 


on body y u + y x2 + vft onset = 

on body C 2 ^ 21 + f 22 + ^ onsQt = 


Whore 1st subscript of denotes the field point location 
and 2nd subscript denotes the location of i.'iPS on circle C-^ 
or Cp causing the stream function. 


We 


also observe that 

^j-'l onset ~ I 1 

^ 2 onset ~ r 2 


sin 

sin 


(3.5) 
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'fll 5 ^'12 5 ^'22 and f> 2 i are calculated as per eqn... (3.3) 
N 

^ = Is. 

N 

^12 r 


Sn(n) 


A 


n=l 


l(n) 


sin n9, 

u n 

sin n0^ 


n 


(3.6) 


bl 


likewise relations of*p 22 and y^ 21 are calculated for field 
point P (r,0) on circle C 2 . By substituting expressions 

for s ^2 an< ^ ^'1 onset ^ rom 3.5» 3.6 in eqn. 3.4 we obtain 
the resulting equation involving 2N unknowns. We therefore 
solve for these 2N unknowns A 2 ( n ) by imposing boundary 

conditions at 2N control points (The resulting linear algebraic 
equation are solved by matrix inversion). 


3.4 Mathematical relations for V gl , V ^ 


V 


rl 


Referring to figure 5(>) 

- JL. P Y11 " HT2 ~ 


1 L ?® 1 


v ©q 


1 1 onset 


•TO © 


1 


(3.7) 


T>vi-12 '~^YL2 


d0. 


:V12 


dr 




bl 


(3.8) 


finally ( 
1 


V 


rl 


n 


nxA 


m(n) 


n 


nxA 


1 n=l 


n 


x cos n©-, + 


1 


1 


n=l 


^ ^ 


n 


bl 


a 
n= 1 


“*l(n) sln n0 l dr bl 

— x 1 + COS e 


bl 



1 


(3.9) 



22 * 


V 


Oil 


~^t > Total _ r Sf‘11 (> 


^ r i 




r-^r- 


+ sin Qn 


n 

H- "* 
n=l r 


nxA 
n+1 


ffl (n) sin n9 l “ mA l( n ) sl ° B h ., -“bl 


or. 


1 


+ X. 

n*l 


"n+T 


Hr 


bl 


1 


n nxA-, r n n cos n0-> d0, 

- ^ ■ ^ (B ll 1 * - sin ®i 


n=l 


bl 


1 

( 3 . 10 ) 


(3.11) 


Similar relations, are obtained for V Q -jr> ang V r2 . 

3*5 Incident flow 

When flow is at an angle of attack o< to the axial direction, 
the resolved component along z direction is [fig. 6b) refers). 

The MPS are aligned along z direction. 0 is measured with reference 
to z direction onset is given by eqn. 3*11. 

^ onset = - r ]_ sin P on Cl 

^onset = T^Ca-r ^ -sin Q^'s P on C2 

W.^ = s in^ = xo(, Taking 1^=1 , max=.l radian, 

Ta^ = . 10 (max . ) 

The relations (3.4), (3.5), (3.9) and (3.10) remain unchanged. 

Since geometry of incident flow is different than that of 

yawing flow, only geometric relations between 0, 0, r , differ. 

d0q ^ r ai 

By substituting proper values of gg- , — ■gjr— in 

relation (3-9) and (3-10) V pl and V 011 are calculate^. 
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3 . 6 Com b ine <d so lution of axiai^ an d lateral flows 

After finding out the perturbation velocity components 
like u x , v r and Vg a from the axial solution and ‘Vg^ and V r 
from the lateral solutions the resultant velocity V n is easily 
calculated by adding the various velocity components in vector 
form. The pressure coefficient Cp (©,x) side force gradients 
dc . dc 

Hx^ 9 Hx^" 9 side force co-efficients C^, C 2 and moment . 
coefficients C , C are calculated as follows. 

Ill v LLL£» 

v n = W* + u x' 1 +■ v r 3 + v © K 
where u x , v r , Vq are axial, radial and tangential velocity 
components ^Ref. fig. 8) of the combined solution such that 

V © = V ©a + V ei 9 v r = v rb ^ V r 


V = 
n 


/ 


T 


~T 


(U + u x ) + v r . 4 . V 


Cp (0,x) ~ 1 - 


V 


n 


Taking =■• 1 


0 2 ..2 2 \ 
= - \2 u v + v„ + \ Q + u v ; 


X 


0 - 12 ) 


dc 

Side force loading 

dc„ 


n f 2 A" 

J- 

S" ’ 


formal force loading 


dx 


i r 

S r ) 


-J 0 

2 A 


r Cp «os & d@ 
Cp sin ©*r d© 


S = ref. cross-sectional area = 2xi$ 


(3.13, 

(3.14) 



Side force co-efficiency c = 
Normal force co-efficient C z = 
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t 

{ do 

j 

(3.15) 

dC 2 x 

j lH3T-^ dx 

(3.16) 


Moment co-efficients (about nose) C 

m Z- 


f 1 

I (+x) (^) dx (3.17) 
■ o 


i do 

C mfe.= ( t' x > ©gf> dx < 3 - 18 > 

. J o 

Analytical solutions of equations 3.13 to 3-18 are difficult, 
so their solutions have been obtained by numerical computations 
by Simpson’s Rule. 


3 »7 Presentatio n of results 

Results of V Q1 and V r in yawing flow over two identical 

and different circles are presented in fig. 21, 22 and 23. 

Results of VM-, and V as function of © in incident flow are 
b± r 

presented in fig. 24. The side load and normal loads of yawing 
and incident flow are presented with the results of axial.' flow 
in fig. 25- Normal load distribution for incident flow are 
presented in fig. 26. 
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Discussion of results and conclusions 


In this chapter the results obtained by the present 
method of representing axially symmetric bodies are reviewed 
and compared 'with those obtained by other methods and by 
other authors. The error in representing the body with and 
without interference has been compared along with distribution 
of pressure and side forces. The factors on which the sensi- 
tivity and accuracy of 3-dimensional dipole method deipends 
namely, calculation of delta-line, choice of form of source 
and dipole distribution and the degree of polynomial are 
discussed? relative significance of lateral flow studied. 
Finally, the suitability of the present method assessed. 

4 • 1 C ompara tive study of results 

As stated in chapter 2, the proposed method of represent- 
ing the axially symmetric bodies by polynomial distribution of 
sources and other dipoles, is applied to the study of flow 
over a family of sharp-ended bodit-s generated by arcs of 
parabola. According to flow tangoncy condition the normal 
velocity should be zero at all points on the surface of the 
body. This condition is exactly satisfied at chosen control 
points but at other points, this condition in general is not 
exactly satisfied. The normal velocity ( v norm ) £ t test 
points chosen midway between the control points called the 
leakage velocity, is a measure of error associated with the 


26 . 


method. ^ norrn normalised with respect to thickness ratio 
has been presented as error in figs. No. 11 to 13 . 

It may be seen that the maximum error i.e. 4.9/. corresponding 
t ° ~&± = • 20 , 7^2= . 2 u and PR = 1.1, 0=180° is comparable with 
the maximum error associated with SDO model of Christopher and 
Shaw. (4.9/0 . The SDoO model of Christopher and Shaw 
gives slightly less error i.e. 4.0/.. 

It is worthwhile to state that all the error results are 
plotted for 0=180° because error at 0=180° is found to be 
maximum. Since boundary conditions, for Obtaining the solution, 
were imposed at 0=0°, where the error was found to be less than 
l/» at all points on checking. The variation in perturbation 
velocity components u , V~ has been presented in fig. 14 to 16 
as function of x, PR, and 0. Results for comparison are 
not available. However, our observations are as follows 

(a) u rises towards - U steeply between x=.05 to x=.00 
because at x=0, there is a stagnation point where total 
velocity is zero i.e. U^+u =0. So u^= -U^. This trend 
is clearly exhibited by the results. 

(b) u appearsto bo proportional to T for small values of XT 
but rises faster for high t . 

(c) Thors is continuous acceleration of flow from nose 

to the ruidele and deceleration from the middle to tail. 
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(d) The variation in ux is maximum along meridian plane 
(9=0) but tnere is no change in u x . clue to interference 
at ©=180 . It shows that interference effects increase 
from 6=180'" to 0=0 from minimum to maximum. 

(e) u x approaches °° in case of touching bodies at x=.50 as 
may be expected. 

Pressure distribution results are plotted in fig. 17 
as function of x, © a't~ER=1;.1 . Cur observations 

ci re as follows • 

(a) Cp variation is maximum near middle porti..n of the body 
along 0=0. The figures show that suction is created and 
attractive loads are acting on interfering bodies in the 
axial flow. 

(b) Interference effects are less near nose of bodies (Cp 
values for isolated and intorferring cases are close) 
because of larger local proximity. 

(c) Cp values at 0=180^ are same for interfering as well 
isolated body. (No interference at 0=180°) . 

(cl) Cp values arc finite near nose. 

(e) Larger loads are created as TJ increases. 

dc v 

Side force loading results for axial, yawing 

and incident flows are presented in fig. 18 to 19 

are compared with the results obtained by panel 


Those results 
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method and those by Christopher and Shaw. It may be seen 

that the results of present study agree with the results obtained 

by’ panel method' and 'Christopher and Shaw'better in close 

proximity near middle portions of the bodies. Our results 

predicted near (x=.lu to .35) slightly over-estimate the 

dc 

side force gradients . The results of normal loading 

Lt X 

for a combined axial and incident flow corresponding to 

a=.10 radian, are plotted in fig. 26. It, may 

be seen that in the range x= . 2 to .35 our results overestimate 
dc 

as compared to the results of panel method and Christopher 
ane. Shaw. But the results are comparable near nose, tail and 
middle portions of the bodies. Results of perturbation 
velocity components V &1 and V r of the lateral flow arc 
presented in fig. 21 to 24. Results of lateral flow 

are not available separately for comparison; however, our 
observati .ns arc- as follows 


(a) Yawing flow produces repulsive forces in the bodies but 
incident flow produces attractive forces as the axial flow 

Cl 2/ 21 $ o 

(b) The flow in general is retarded in yawing case whereas the 
incident flow is accelerated. 

(c) V (V ) as stated earlier is a measure of error, the 
' r norm 

values .,f V ■ are large in incident and yawing flew. 
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At PR=1 . 1 , identical bodies 
V = 2U x 10 ^ incident flow 

XTii 

—2 

Vnorm = 15 x 10 yawing flew 

These values are with reference to onset lateral flow 
i.e. VV^ and V^. Even though, the errors are large with 
reference to lateral flow but these errors are small and of the 
same order as the axial flew when compared with i.e. 

m =2xl0 , 1.5x10 for incident and yawing flow as 
compared t- U^. The contribution of lateral flow is very 
small because angles ,f incidence and yaw are small (limited to 
.10 radian) as dictated by slender body theory. Side force 
(^-.-efficient (Cyj -detained by integration of the gradient 
-■rf“ results are presented as function of PR and Xj in 

C.i X 

fig. 20. 

This suggests that 

One con make a choice between the load encountered at a 
particular proximity if external stores in those proximities 
are to bo carried. It is seen that the side frrcc- co-efficient 
steadily increases with proximity and ‘ . > until 

PR=1.5. The value of rises steeply between PR=1.5 tc 1. 

4 « 2 Sensitivity, improv eme nt and application of the method 

The accuracy _-f the method, depends on accurate determina- 
tion of delta line by taking 30-approach, choice of polynomial 
functi.ns f (x) , fd(xj and the degree of polynomial. In the 
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present study delta line calculations are based -n 2D 
approach. This node’s improvement for greater accuracy. 

Regarding proper choice of f(x) and f^(x) we have taken into 
consideration the experience < f other authors and consequently 
w c could represent body (isolated case and with interference) 
by a polynomial of degrees 8 for the profile chosen. It 
must be mentioned here that a polynomial of higher degrees 
than 8 did not improve accuracy; rather polynomials of degrees 
greater than 10 led t^ ill-conditioning - of matrices. 

This method is applicable t:. sharp nosed bodies but it 

can be tried v-ith blunt nosed bodies by redistribution of 

polynomial functions between x=.02 to ,98 or so. This method 

has been tried on identical, symmetric, parabolic bodies with 
T. 

varying '"between 0,10 to 0.40 ( =?10, 20, 40/,) and for different 

proximities . The body profile being represented by r(x)=2Z3 
( x-x^) , where 17 =thickness ratio (d/l). The method gives 
good results for interference between bodies of unequal 
thickness ratios (17) also. The method is easily and readily 
extendable t.. non-identical bodies laterally separated, 
staggered and with modifications to skewed bodies t^o. Our 
further suggestions for extension and application of the 
method arc as fellows. 

1. Interference between axi-symmetric bodies in compressible 
subsonic and supersonic flows. 

2. Interference between cluster of axi-symmetric bodies. 


4.3 


(a) 


:ENTR/V. 1 -SR ARY 


Conclusi- ;_ns 
Following exclusions 
j^ isc rete^ distribution 



arc drawn from the- present 
of singularities 
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study . 


Representing an axi-syrnmetric body by discrete 
distributions as largo as 35 failed because of large 
leakage velocities at test points. This method was given 
up after initial trials. 


(b) This study suggests a fairly accurate method of rapid 
do termination of interference effects between axially 
symmetric bodies. This may supplement more accurate 
panel methods by providing an approximate initial 
solution over which panel methods may be efficiently used 


improve accuracy. 
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Evaluation of integrals / "- ^ 7 


-372 


I m = * 


-m 

x dx 


x Tl ^ , ( 2m- ' 3 ) x ^_ jx m "^dx 


R' 


V 2 = 


R' 


s/r 


(m-l) (x 2 +r _/ 
(m-2) 




J 


x m " 2 dx 


R 


37T 


— 9 o 

At m=2 where R=(x-x) +r 


-x 

Tk 


+ x S X ,3j2 dx + / 


tf 


■§--“+ x x- I 1 + I Q 




R' 


t _ r M,, _ , rsh tx-x, 

o ~ J \jh ~ rtrSu 


I, =/-^ 


R 
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Selutien ef the _co-eff icients ef the polvnemial distribution 
functions 

From flow tangency condition we have 


v 


s2 


Aj9 

r 

dr(x) 


■^u 


dx 



X 



e A v r = 

V s2 

+ v di - 

v a2 

^ u x = 

U s2 

+ u di + 

u d2 

and u ~ 
s2 

are 

the induced 


(2.15) 


presence of body in the flow represented by a function f^Cx) . 
Function f^x) is calculated similar to the method given in 
para 2.3 (eqn. 2 . 11 ) ■> So v g2 and u g2 are known parameterw. 

Subs titutino for v and u w in eqn (2.15) we obtain a 


relation 


dr( x; 


^ u dl +u d2^ ox "\ v dl 


"d2 


y 


s2 


s2 


dr( x ) 
dx 


(2.26) 


» v dJ , v^^ i n terms of f 2 d^ ate. as 


1 f ld (x) <J5 


R 


Putting u dl , u d? 

given in eqn. 2.17 and 2.18 

r 1 f ld (; ) dL r 
^dl - r dll { - q< 3/2 “ : 

where R^.= (x-x)^ + r dll^ Fef . fi<3° 4 

n 

2 

j=o 


3/2 


(2.27) 


n 

v,, = \ r . , i x 2 

d 1 7 ci 1 1 j _ o 


1 '(jj - r x x S °d(j) 


0( 


•dtj) 


A 

J 


(2.28) 
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By moving the field point P at (n+l) control points 
we obcain (n+l) eqns , and 2.28 is modified as follows 
in vector form. 

v dl^ = -) r dll^^ x x cos (© b -G) -- Aa ( i , j ) x r- L (i)/ x c o(j) 

'' J (2.29) 

By resolving velocities in radial direction along ©. 


Similarly 

v d2(i) = _) r dl2(i) x B d(i,j) x cos r bl(i) xcos(0,--e)l 

l . 

where Acl(i,j) : 


n 

£ 

3=o 


'ci( j) 


'b( j) 


J 


x Cd(j) 

.1 553 H v (2-3°) 


o li 


35 J Cix 

'3/2 


3(i» j) 


n 

£ 


j=o bd(j) 


I bd( j ) ~ ^ 


" x J dx 


o R 


bd 


3/2 


Bg(i, j) 


n 

£ 

j=o 


'b( j) 


j - ^ 

Further substitution gives us 
vd x £i) = Wj^Cijj). x cd( j) 


* R b ~ ^ x ” x i^ 2 +r bl (i) 2 


R bd = (x “ x i )2 + r dl2( i) 2 


(2.31) 


Vd 2 (i) = VV 2 (i,j) x cd( j) 
where W^i.j) = r c m(i) 
W 2 (iJ) = r dl2(i) 


X 


Ad (i ? j) x cos (©8-©) - AA(i,j) x r^i) 
Bd( i , j ) x cos (^-G) -8$(i, j) x rb x (i)x 

X COS'^ - ©/ 


X 
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Similarly ucL (i) = J 


1 f ld (x) ( X j~ X jdx r 1 f ld (x)(x i" x) dX 
1.3/ 9 ' “ */ ““ 9/0 


o r; 


= x i 2_ Id(j) -Z Id( j+l) -/x i Z I. - Z l(j+l)j 


j=o y 

(2.32) 


which in Matrix Form 


( i) _ E -j xAD j i ,, j ) — Ad(i, j+l)^— jx^xAA(i,j) — AA(ijj+l)^ J 

l 


xcd( j) 


ud 2 (i) = [ ^ c i xBd(i , j) - bd(i, j+l)? - ^p^xBBU, j) - BB(i,j+l)X ] 

xcd(j) (2.33) 


on substitution, 


ud 1 (i) = UU i (i,j) x cd(j) 
ud 2 (i) = UU 2 (i,j) x cd( j) 


(2.34) 


where UU,(i,j) = jx- x Ad(i,j) - Ad(i,j+l)/- ix- xAa( i , j)— AA( i » j+l) ( 

1. J I 

UU (i,j) =^x. x Bd(ijj) - Bd(i, j+l) L- jx. BB(i , j)~BJB( i, j+l) l 


J L' 


(2.35) 


hqn. (2.41) can be put in matrix form as 

dr. 

C gJU 1 C±, J) + UU 2 (i,j) (x -|w x (i,j) -VV 2 (i,j)j] X cd(j) 


r dx 


MHlVC ^ Ct> ‘ v =f' "tfu J" c \ 


= v. 2 (i) - u g2 (i) x 


(2.36) 


FX = B 


X = F B 

where B = 8(i) = K.H.S. of (2.36) 7 
X = X(j) = Cd(j) unknown. 


(2.37) 


F==F(i,j) = 'UU 1 (i,j) + UUgli.jHx -WgU.oJY 

? J L ( 2 .i) 
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Appendix 1 C* 
Geometric Relations in 


Flow 


Yawing Flow : (Ref. fig. 5b) 

The main relation from which other derivatives and 
eqns. are drawn. 

r^ sin 0 1 = R fe ^ x sin (it - 0^) 
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a x x sin G^ 
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(Ref. fig. 6b) 
a r, cos ©. 
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o Source C = -20 Body eqn. 
A Dipole C= '20 



Fig. 10 Strength distbn. of source & dipole PR=M 
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Fig.12 Error variation with proximity 
1= -20,0 = 180 
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Fig.13 Error variation with thickness 
PR = 1-1,0 =180° 
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over surface 
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Fig.15 Variation of axial pert bn. velocity with 
' proximity Z)-Zf20 , 0 = 0° 
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Fig. 16 Variation of tangential velocity over surface 
Ci=C 2 = -20 , PR = H 





Panel - method (35x20) 



Fig-18 Side load distribution with x in axial flow 
Cr L 2 = -20, PR =11 





Results on body C-j 



g.21 Variation of tangential velocity with 9 
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F i g- 22 Error v q r i g t i o n with proximity 
Vnorm = Vr = Error 




° VQ{ , Interference results 
A Vgf , Isolated body results 
O Vnnrml^f fOf) 





A Axial flow 



Fig.25 Axial and lateral flow combined 
Ci = C? =-20, PR = M 



